Mass gap in scalar glueball model 
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On the basis of a nonperturbative scalar model of gluon condensate the model of glueball is 
considered. Two scalar fields describe quantum fluctuations of gauge potential components belonging 
to a small subgroup SU(2) C SU(3) and a coset SU(3)/SU(2) correspondingly. In this consideration 
nonperturbative quantum corrections appear. The corrections are connected with a hidden structure 
of operators of strongly interacting gauge field. A regular solution for two scalar fields is obtained. 
The solution is a ball filled with quantum fluctuations of SU(3) gauge field. The space distribution 
^J ■ of the dispersion of quantum fluctuations of gauge potential is presented. 
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Non-abclian gauge theories are central to our current understanding of physical phenomena (excluding gravity). 

j~h, The perturbative analysis of such theories is fairly well understood by now, having been extensively developed over 
the last three decades. Many of the nonperturbative aspects are also more or less understood at a qualitative level. 
However, it is fair to say that, as of now, we do not have calculational techniques or detailed understanding regarding 
nonperturbative phenomena in non-Abelian gauge theories. 

Let us summarize the standard lore attributed to quantum chromodynamics (QCD): (a) QCD must have a mass 
gap, i.e., every excitation above the vacuum state must have energy A > 0; (b) QCD must exhibit quark confinement, 
i.e., the observed (physical) particle spectrum are color SU(3) invariant, despite the fact that QCD is described by an 
underlying Lagrangian of quarks and non-abelian gluons, which transform non-trivially under color SU(3) symmetry; 
fT} ' (c) QCD must exhibit chiral symmetry breakdown, i.e. the vacuum is only invariant under a subgroup of the full 

symmetry group acting on the quark fields, in the limit of vanishing current quark masses. 
f— ^ ■ Item (a) is essential to understand the short-range of the nuclear force. Item (b) is essential to account for the 

absence (unobscrvability) of individual free quarks. And finally, item (c) is essential to justify the spectacular current 
algebra predictions. 

The perturbative description of elementary particles is essentially based on the field-particle duality which means 
that each field in a quantum field theory is associated with a physical particle. On the other hand, the situation in 
/\ ' QCD is more complicated, however. For a description of confinement of quarks and gluons within the framework of 
H \ local quantum field theory, the elementary fields have to be divorced completely from a particle interpretation [l[ . 

Therefore, to study the infrared behavior of QCD amplitudes non-perturbative methods are required. One promising 
approach to non-perturbative phenomena in QCD is provided by studies of truncated systems of its Dyson - Schwinger 
equations, the equations of motion for QCD Greens functions. The underlying conjecture to justify such a truncation 
of the originally infinite set of Dyson - Schwinger equations is that a successive inclusion of higher n-point functions 
in self-consistent calculations will not result in dramatic changes to previously obtained lower n-point functions. To 
achieve this it is important to incorporate as much independent information as possible in constructing those n-point 
functions which close the system. 

Usually the attempts to solve problems in QCD are based on notion of gluon - a quantum of SU(3) gauge field. 
Tacitly assuming that quantized field is a cloud of quanta (gluons in QCD). But this is not evident for a quantum field 
theory with strong interaction where it is necessary to use a nonperturbative quantum field theory. It is possible that 
strongly interacting quantum fields more similar to a turbulent fluid where is a preferred direction (for example, it can 
be flux tube stretched between quark and antiquark) and fluctuations around nonzero field there exist. In standard 
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quantum field theory the interaction between quanta happens in a vertex. The vertex is O-dimensional object where 
O-dimensional objects - particles do interact. All it gives rise to the idea that strongly interacting quantum fields can 
not be described on the language of particles - quanta. Mathematically it means that in this case it is not possible 
to use Feynman diagram technique to the description of strongly interacting quantum fields. The consequence of 
this statement is: if we can not use Feynman diagram technique then we can not use the notion of quanta for the 
description of strongly interacting quantum fields. In other words strongly interacting quantized fields do exist but 
they do not formed with quanta. For example, quantum gluon field does exist but gluons (in some energy region) do 
not exist, for gravity this statement is even more strict: quantum gravity does exist but gravitons do not exist. 

Here we use the idea that the operator of strongly interacting field can be decomposed on constituents (hidden 
structure of the operator). In order to these constituents are invisible we assume that they are nonassociative operators. 
By calculation of expectation value we need to rearrange of brackets that leads to the appearance of additional terms 
in the same way as Planck constant appears by the operators permutation in standard quantum field theory. Following 
this way we calculate 4-point Green's function. In this case additional terms appear in the consequence of the operator 
hidden structure. In order to avoid solving an infinite set of nonperturbative equations connecting all Green's functions 
we calculate gluon condensate assuming that: (a) 4-point Green's function approximately can be presented as a bilinear 
combination of 2-point Green's functions; (b) 2-point Green's function can be approximately presented as the product 
of scalar fields (scalar approximation). It means that these scalar fields describe the correlation of quantized fields in 
two separated points. If two points coincide then the scalar fields describe the mean value of the square of quantum 
fields, i.e. the dispersion of quantum fields around their vacuum value. In order to obtain equations for the scalar fields 
we average the SU(3) Lagrangian (which is none other than gluon condensate). Finally we obtain the Lagrangian for 
the scalar fields. We solve these equations numerically in the spherically symmetric case. As result we obtain a ball 
filled with fluctuating quantum SU(3) gauge field that can be presented as a model of glueball. The solution does exist 
for a special choice of the value of the scalar fields at the center only. One can say that corresponding scalar fields are 
eigenfunctions of corresponding nonlinear differential equations. It means that the ball filled with fluctuating SU(3) 
gauge field can have unique value of energy. The ball can not have the energy value small as possible, i.e. we have a 
mass gap in this model of glueball. 

Thus the main idea considered here is to offer an approach to the quantization of strongly interacting fields do not 
using the notion of quantum. 

Some reviews about nonperturbative physics in QCD arc: (a) the studies of QCD Greens functions and their 
applications in hadronic physics in (l|; (b) light-cone quantization of gauge theories in [2J; (c) the theory and phe- 
nomenology of instantons in QCD in J3|; (d) numerous aspects and mechanisms of color confinement in QCD arc 
surveyed in [4| ; (e) in Ref. [5| recent results in the physics of glueballs with the aim set on phenomenology is reviewed 
and the possibility of finding them in conventional hadronic experiments and in the Quark Gluon Plasma is discussed. 
The nonperturbative approach for the quantization of strongly interacting spinor field was presented in Ref. [6( . In 
Ref. [3] general arguments are presented that only scalar glueballs should exist using the requirement that the gluon 
field have to satisfy the canonical angular momentum relationship. 

II. SCALAR MODEL OF GLUEBALL 

Following to calculations from Appendix [A] we can obtain an effective Lagrangian describing non-perturbative 
quantum fluctuations of SU(3) gauge field. The Lagrangian is obtained by averaging of SU(3) Lagrangian using 
a hidden structure of field operators. The hidden structure means that the operators of strongly interacting fields 
have can be decomposed into the product of some operators (constituents). The constituents should be invisible 
(unobservable) . Such invisibility can be only if the constituents are non-associative operators. The effective Lagrangian 
describing gluon condensate is 

Ceff = \ (d, X f + \ (d^f -±(X 2 - m\) 2 ~^{<P 2 - ml) 2 - ^V + const (1) 

where the scalar field \ an d 4> describe quantum fluctations of gauge potential belonging to a small subgroup SU(2) C 
SU(3) and coset SU(3)/SU(2) correspondingly. It is necessary to note that for the construction of a scalar model of 
glueball these fluctuations arc essentially different. 

Our goal in this section is to find spherically symmetric solution describing non-perturbative quantum fluctuations 
of the SU(3) gauge field. The field equations for the Lagrangian ([T]) are 

□0 = -0 [A 3X 2 + A 2 (0 2 - mj)] , (2) 

□X = -X [A 3 2 + Xi (x 2 - ml)] . (3) 



Since we are searching for a sphericlly symmetric solution then the equations have following form 
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0" + -0' = [A 3X 2 + A 2 (0 2 - ml)] , 



x" + -X' = X [A 3 2 + Ai (x 2 - ml)] 



(4) 
(5) 



where prime (•••)' = d(- ■ ■ )/dr. The numerical analysis shows that regular solutions exist for a special choice of 
boundary conditions 0(0), x(0) only, when the parameters Ai ; 2 and mi^ are fixed. One remark is that technically 
(for the numerical calculations) it is easier to fix Ai,2 , <^(0), x(0) and search special values for m\fi. For the details 
of numerical calculations one can see Ref. Q • 

By some technical reasons we will fix 0(0), x(0) and using shooting method we will find special values of mi. 2 for 
which regular solution does exist. We will search regular solutions with the following boundary conditions 
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After that we redefine functions and parameters in Eq's ((4]) ([5]) in following way: 
rfi2 = raijra\^ x = m\r and fh\ = m\jm\ = 1. Then we have following equations 



0" + -0' = 
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x" + -x! - 
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= X 



A3X 2 



A 3 2 + Ai (x 2 



1 



(0) 

(7) 

(r)/mi, x(x) = xW/mi, 



(8) 
(9) 



here (•••)' = d/dx. Such redefinition allows us to consider mi. 2 as fixed parameters and equations ([5J (j9)) as nonlinear 
eigenvalue problem for eigenvalues 0o,Xo and eigenfunctions 0, \. The results of numerical calculations are in Fig's 
[T] and [5] presented. For the numerical calculations we use following values of parameters Ai = 1.0 , A2 = 0.1 and 
A 3 = 1.0. 

The asymptotical behavior of the functions 0(x),x(x) is 
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TABLE I: The values of 7711,2, E as the functions of 0o, Xo- 
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TABLE II: The values of 0o, Xo, E as the functions of rhi- 



The gluebal mass m g can be calculated in following way 



mgC - =E = 4nJr*e (0, X ) dr = 4™* J x*e ($, X ) dx = mi E. 





(12) 
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FIG. 1: The curves la, 2a, 3a and 4a are the profiles of 
x(x); lb, 2b, 3b and 4b are the profiles of <f>(x) for mi, 2 
given from the Table [TT1 
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FIG. 2: The curves 1, 2, 3 and 4 are the profiles of di- 
mensionless energy e(x) for mi,2 given from the Table HT1 



The energy density is 



e = \<P' 2 + x' 2 + Vtf, X) = m\e = m\ (\~^ + X ' 2 + V$, x) 



(13) 



In Fig. [3] the values of (j>o,Xo,E from the table [TTI and corresponding fitting curves (fTl| - (fl"l))) are presented. 
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FIG. 3: The values 0o, Xo, -E 1 from the table ITT1 and fitting 
functions from Eq's (fT4"l) - (fl"6)) . 



FIG. 4: The profiles of the dispersions T> (A M ) and 
V(A™). 



The fitting functions of <j)o(m2), Xo(w 2 ) and -E(m 2 ) are 



</>o(m 2 ) = 1.001 -0.4689m 2 + 0.05063m2, 
Xo(m 2 ) = 0.000612 + 0.9804m 2 -0.40428m2, 
E(m 2 ) = -0.0000967 + 0.01356m 2 + 0.02146™^ 



(14) 
(15) 
(16) 



The numerical analysis shows that 

m^ 2 — > 0o = 1, (17) 

E -^» 0. (18) 

It is interesting to see how the dispersion of quantum fluctuations 

V (A?) = G 2 (x,x) = {A™(x)A m »(x)) ex m| - 2 (x), (19) 

V (4) = G 2 (x, x) = (Al(x)A»»(x)) ex x 2 (x) (20) 

in the scalar model of glueball in the space are distributed. These profiles in Fig. |4] are presented. The negative 

value of T> (A™) means that the vector A™ is space-like one (we use following Minkowski signature: (+, — , — , — )). 
The negativity and positivity of the dispersions V (A™) and V (A h ) means that in the first case the components 
A™, 1 = 1,2,3 are dominant and in the second case the component A b t is dominant. 



III. DISCUSSION AND CONCLUSIONS 

The main result of this research is that we have obtained an approximate model of gluon condensate using two 
scalar fields approximation and the spherical symmetric distribution of gluon condensate as glueball is interpreted. 
The numerical calculations show that for given parameters ih\ 21^1,2 there exists a unique value of scalar fields at 
the center 0(0) and x(0) that the solution is regular. It means that in the theory with given parameters mi,2,Ai.2 
there exists a unique excited state of nonperturbativc quantized SU(3) gauge field with the energy E{rh\ y 2, Ai^)- This 
unique value of the energy is a mass gap A (mi,2, Ai.,2) = E (mi,2, Ai^) for given parameters 7711,2, Ai,2- 

Let us remember that scalar fields \, describe quantum fluctuations of the components belonging to small subgroup 
and coset. Thus our solution describes a ball filled with quantum fluctuations of the SU(3) gauge potential. The 
quantum fluctuations of A\ e SU(2) and A™ € SU{3)/SU(2) have different behavior. 

Interesting problem in two scalar fields approach is its gauge invariance. The condensates (|A14|) and (|A15[) are not 
gauge invariant. The question is: is it necessary to this approach or it is possible to avoid this problem ? We see two 
possibilities: the first one is that the nonperturbative quantization leads with the necessity to the violation of gauge 
invariance of gauge theories; the second one is that one can use a minimal value of the condensates (|A14|) and (|A15|) 
following to Ref's [§] [13]. 

Let us note that a discrete spectrum of the energy of localized field distribution have been found in Ref. [LjJ for 
non-linear Klein Gordon and Dirac equations. The difference with our case is that non-linear Klein Gordon and 
Dirac equations have infinite energy spectrum where the first state can be identified with mass gap. In our case two 
interacting Ginzburg - Landau equations © © are constructed in such a manner that we have only one regular 
solution which we indentify with the mass gap. 

The most unexpected feature of presented approach is the appearance of a hidden structure of operators of strongly 
interacting fields. In order to mask this structure it is absolutely necessary make the constituents of the operator 
decomposition as non-associative quantities. Such non-associative decomposition leads to nonperturbativc quantum 
corrections which appears as exponential decreasing of quantum fluctuations, see (|10p (jTTJ) . The values of these 
corrections can be in principle estimated as the size and energy of glueball. 

Another interesting feature is that the components of gauge potential belonging to a small subgroup SU(2) c SU(3) 
and coset SU(3)/SU(2) behave variously: in quantum fluctuations of gauge potential belonging to the SU{2) the A t 
component is dominant but for quantum fluctuations of gauge potential belonging to the coset SU (3) / SU (2) the A 
space components are dominant. 

It should be noted that the parameters 7771,2 in Eq's ((2|) <j3j> can be omitted and obtained in following way: (1) scalar 
fields x an d <t> are complex-valued fields; (2) the parameters 7771,2 in Eq's © $5§ are removed; (3) the solution of Eq's 
(@| (0) (without 777i. 2) is searched in following form 4>{t,r) = e _I " 2 *0(r) and x(^ r ) — e~ lUlt x(r). After that we have 
the same equations as (QJ © but with the replacement mi, 2 — > ^1.2- It means that mathematically both approaches 
are equivalent. But physically these approaches are not equivalent. If we work with u>i .2 then the parameters of the 
theory are Ai,2 only, and wi,2 are eigenvalues, i.e. for every 0(0) — > 1 and x(0) — > there exists the regular solution 
and the glueball mass decreases to zero. Consequently the mass gap does not exist in this approach. But, as we have 
shown above, for mi, 2 approach the mass gap does exist. 
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Appendix A: Nonassociative operator decomposition — hidden structure of field operators 

In this section we briefly describe nonassociative decomposition of strongly interacting quantum field operators (for 
details see jl3|). 

Let us assume that operators of SU(3) gauge field fields A B (x u ) can be decomposed into non-associative (n/a) 
constituents e), and $ ?s : 

K = 4^ B , (ai) 

here i is the summation index, /i is spacetime index and B == 1,2, ••• ,8 is the SU(3) color index. Although the 
constituent operators e* , $> lB are not associative, the basic idea requires their product to model an associative operator. 
In more mathematical terms, the operators e % and $ lB are elements in a n/a algebra A, i.e., e l , $ lB <G A/G, which 
contains an associative subalgebra G C A, such that A B = e* $ lB € G. It is necessary to note that the operator A B 

models observable quantities, whereas the n/a e! I ,$* B are unobservable. In classical case the decomposition similar 
to (|A1[) is applied for SU(2) spin-charge decomposition [l3[ - [Tg|. In solid matter the decomposition similar to (|A1[) 



is applied to fermion operator and is known as slave - boson decomposition [17| - [19[ . In both cases the constituents 
are associative ones. 

The SU(3) Lagrangian is 

C = F B * V F* V (A2) 

where F B V = d^A 3 - d v A B + gf BCD A^A^ is the field strength; B, C, D = 1, . . . , 8 are the SU(3) color indices; g is 
the coupling constant; f BCD are the structure constants for the SU(3) gauge group. 

Our goal is to average quantum version of the Lagrangian (|A2|) which is none other that gluon condensate 



<£> = {F B ^F B V ) (A3) 

with some approximations concerning to 4-point Green's function 

G B ^ p DE (x 1 ,x 2 ,xs,x 4 ) = (A B ( Xl )AC(x 2 )A D (x 3 )A E (x 4 )) (A4) 

where (■ ■ • ) is a quantum averaging over some quantum state \ip), and A B F B in (|A3[) (|A4[) arc the operators. 

In order to calculate an expectation value we have to define the action of operators on a quantum state. We will 
require: 

a b m = (4$ ?B ) iv) = 4 (V s m) (as) 

this is the same rule as for the associative case. But when having two or more nonassociative operators, there is 

((4 $?fl ) 4) {& c w) + (^C iv>> = (A6) 

where the associator (rh 2 ) is defined as follows 

(e^ B ) (el& c ) = (((4.^ s ) eP^'A + (m 2 f^ . (A7) 



,BC ■ 



One can say that (m2) is the second Planck constant but with different dimension. The same can be done for 

A 3 m 



A B A?A? |V.) = eU & B (ei(& c (4 (**" |^>)))J J + {mrf? A* \1>) + (fhsf™ |V>) 
here we use following non-associativity 

For A 4 \tp) 
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A B A C V A°A E |^> =< $ iB ei <^ c e* 4> feD (C (^ |^))) 



(A10) 
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+ {*»)£ K A » M + (^C K M + (m 4 )^r ivo 



here we use following non-associativity 

(ej,fc iB ) U U' c (ej(^(e* (V 7i£ |V») 
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The next assumption is that the SU(3) gauge potentials A B <G SU(3), B = 1, • • • ,8 can be decomposed on two parts 
having different quantum behaviour. Physically it means that quantum fluctuations of A a <G SU(2) c SU(3),a = 
1,2,3 and A™ <G SU (3) / SU (2) , m = 4, • • • ,8 are different. The gauge potential A a belong to the small subgroup 
SU(2) C SU(3) and gauge potential A™ belong to a cosct SU(3)/SU(2). 
We assume that 
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e^xx) (#*(*!) (ei(x 2 ) (V 6 (x 2 ) 
J,(xi) (^ ro (x x ) (4(x 2 ) (V"(x 2 ) 
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(A12) 
(A13) 



here Cf ll ; mn are some constants and x(x),</>(x) are scalar functions that describe quantum fluctuations of gauge 
potential components belonging to a small subgroup and a coset correspondingly. Thus 

(4(x!)^(xa)) = CfrxfriMx*), 

(A";( Xi )A:(x 2 )) = C->(x 1 )^(xa) + (m 2 )™ r 

here we assume that (jh 2 ) = and consequently 

(4(x)^(x)) = C b Xx 2 {x)^ ai x 2 {x) 
(A™(x)A m »(x)) = C mm ^ 2 {x) + {m 2 ) mm ^a 2 ^{x) + {rh 2 ) 

where cm = C^ and a 2 = C TOm £. Now we can calculate Ud^ 7 " 
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For the calculations for A 4 let us write 

rBCD aC aD rBMN aM[i aNu rabc rade \h \c Adfi Aeu _i_ cy rabc raran \b \c Amp, Any i 

^rrrna nrpq j^rn j^a j^p^. ^qu , rBmn rBpq j^m j^n, j^p^. j^qv 



(A18) 
(A19) 

(A20) 



here b,c,d,e,= 1,2,3 and m,n,p,q = 4,5,6,7,8. Now we need some assumptions about 4-point Green's function. 
Schematically our assumption is that (^4 4 ) ~ (A 2 ) + (^4 2 )- Mathematically it looks as 



4 ( & h ( el[ $ Jc (e^(V d (C (^ e IV)))^ y 



'e*(* a {4 (^ IV-))) )) (e^($ Sd (e™ (<^ e |V)) )) 
The same for the coset components 



si(*>"(e*(#^(e*(** h «|^))) j _ 



ej, ( $ ir " 



(4 (V" |V))) )) (e*(***(e™ (**« |V)) ) 



For the term A h A c v A miL A nv we assume that quantum fluctuations between gauge potential belonging to the small 
subgroup and coset do not correlate 

(A^AIA" 1 " A nv ) = (A^Al) (A^A™) . (A23) 

Thus using ([ATO]) . (JA"P2"T) . (|H4| and ([MT|) we obtain 

(V I / Qbc / Qde A^ A d ^ A e,y | V) = — X 4 (A24) 

where % = / a6c / ode C h ^ v C de ^ v and we use u sual a ss umpt io n tha t (A^ ) = . 
The same is for coset gauge components f (|A10j) . (JA13I) . (|A15|) and JA22p ) 

(V |/ mnp / m9 M"^A^A ri/ | V) = ^ 4 - ym^ 2 + Ci (A25) 



npqr^iu 

For (^ J 4£ J 4 m 'M r " y ) we have 



fmnpfmqr ^ 



(V |/ a6c / om M*^A m M w '| V) = yX 2 2 - ymiX 2 (A26) 

where 4f = f^cjamnQbc^Qmn^^ -km? = jabc jamn (^^"F qbc 

For 

(i^™i w i' 1 '') « (A b ) (A^M 9 ") = (A27) 

as we assume that as usu al (A 2 n +1 ) — 0. 

Collecting all together (|A18|1 . (|A19|) and (|A24[) - (|A26[) we obtain an effective Lagrangian for two scalar functions 
X and 4> describing fluctuations of gauge potential components from a small subgroup SU(2) C SU(3) and coset 

SU(3)/SU(2) 

Ceff = l^xf + lid^f-V^x), (A28) 

Vfax) = Y(x 2 -™ 2 ) 2 + y(0 2 -^) 2 + ^0V + C (A29) 

here we have redefined the scalar fields <f> and x i n such a way that to get rid of constants a\^ in front of ( d^cf) ) and 
(9fiX) i constants Ai.2,wi,2 a re redefined on corresponding manner and constant C is unessential for field equations. 
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